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Preface 


This  translation  of  Chapter  7  of  Ya./Tj,  Al' pert's  book 
"Rasprostvanenie  Radiovoln  i  lonosfer^"  was  supported  by  the 
Advanced  Research  Projects  Agency./ This  chapter  deals  primarily 
with  the  propagation  of  very  low  frequency  radio  waves.  Also  in¬ 
cluded  are  Appendices  6  and  7  of  the  book  since  they  also  deal  with 
long  waves.  Just  before  completing  these  translations  Dr.  Alex 
Nennsberg  passed  away.  Mr.  J.  B.  Reubens  and  I  have  attempted  to 
carry  out  the  editing  of  Dr.  Nennsberg' s  literal  translation.  At  the 
same  time,  I  have  taken  the  liberty  to  add  a  number  of  editorial  foot¬ 
notes  which  might  be  of  interest  to  the  reader.  Mrs.  Eileen  Brackett 
has  prov  ided  valuable  help  in  the  typing  and  preparation  of  this  report. 


James  R.  Wait 
August,  1961 
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Long  Waves 


As  previously  pointed  out,  long  radio  waves  cover  the  range 
between  2,  000  to  20,  000  meters.  HoweveT7  the  fundamental 
peculiarities  of  long  waves,  which  separate  them  from  medium 
waves,  appear  in  practice,  starting  with  waves  between  5  and  6  km 
in  length  and  greater,  at  relatively  large  distances  from  the  source. 

The  height  at  which  the  ionosphere  begins  at  these  low  fre¬ 
quencies  is  composed  of  only  several  wavelengths.  The  transition 
region  at  the  base  of  the  ionosphere,  where  electron  density  and  con¬ 
ductivity  increase  rapidly,  is  of  the  order  of  one  wavelength  or  much 
smaller.  Fur  ihi^  reason,  the  D-region  of  the  ionosphere  is  not  an 
absorbing  layer  for  long  waves,  as  it  is  for  the  rest  of  the  waves  in 
the  radio  frequency  range,  but  serves  as  a  conductive  wall  of  the 
waveguide,  directing  ^’->6  flow  of  radiation,  and,  by  this  effect,  con¬ 
tributing  to  long  distance  propagation.  Long  waves  in  the  earth- 
ionosphere  waveguide  are  excited  as  a  spectrum  of  waves.  Even 
near  the  radiator,  the  interfering  character  of  the  field  structure, 
depending  on  distance,  caused  by  the  superposition  of  these  waves, 
may  be  observed.  At  low  radio  frequencies,  at  a  point  sufficiently 
remote  from  the  source,  direct  and  reflected  waves  from  the 
ionosphere  seldom  appear  separately.  For  example,  at  distances 
above  2000-3000  km,  only  one  sufficiently  intense  wave  remains. 
However,  this  so-called  zero  mode  of  the  spectrum  of  excited  waves 
is  of  a  purely  waveguiding  nature  and  is  not  a  direct  wave.  Hence, 
the  classical  approach,  based  on  the  strict  solution  of  the  diffract- 
ional  problem  (-see  para.  22),  as  well  as  ray-tracing  interpretation 
of  the  picture  of  propagation,  taking  into  account  the  influence  of  the 
ionosphere,  are  often  not  suited  fur  describing  the  propagation  of 
long  waves. 

In  recent  years,  the  theoretical  and  experimental  investigation 
of  long  waves  has  been  considerably  broadened.  A  special  interest 
has  been  shown  in  the  propagation  of  electromagnetic  waves  of  low 
and  super-low  frequency.  In  particular,  such  problems  arise  in  the 
study  of  the  so-called  atmospherics  solitary  signals,  excited  by 
lightning  discharges  propagating  along  the  earth,  and  also  of  other 
signals  of  similar  type.  These  considerations  are  also  important 
in  the  study  of  the  so-called  "whistling”  atmospherics  in  the  same 
frequency  range.  These  are  excited  by  lightning  discharges  which 
leak  through  the  whole  thickness  of  the  lower  part  of  the  ionosphere 
and  then  propagate  along  the  force  lines  of  the  magnetic  field  of  the 
earth  in  outer  space,  traveling  many  thousands  of  kilometers  before 
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their  return  to  the  earth. 


After  an  interval  of  almost  40  years,  interest  in  the  study  of 
long  waves  hWs  revived.  It  is  precisely  during  this  latter  period, 
however,  thara  physical  theory  of  propagation  of  long  waves 

developeJ^which  provides  methods  of  comput^ion  and  makes 
it  possible  to  compare  theory  with  experiment.  tii«ei 

new  experimental  results  obtained  from  tau^sources: 

direct  experiments  with  radio  waves,  and  indirectly  by  using 
atmospherics.  Therefore,  at  present,  it  is  not  necessary  to  be 
limited  to  radio  waves  created  by  radio  arrangements  only,  but  it 
becomes  possible  to  broaden  the  investigation  down  to  elMtromagnetlc 
waves  at  the  lowest  observed  frequency  (i.  e.  ,  10-12  kc/vf. 


Iw  thw  ^ome  results  of  theoretical  computations  are 

tVr  including  super-low  frequencies  down  to  some  tens  of 
cycles  per  second.  The  results  of  experimental  investigations  of 
atmospherics,  particularly  the  whistling  atmospherics,  whisk  arf  nf 
sets iLUt'^rrtiwwtlie"  amps  nf  thift  are  briefly  analyzed.  It 

should  be  noted^at  the  question  of  electromagnetic  waves,  in  the 
range  3-5  kc/rand  lower,  remains  uninvestigated  with  adequate 
precision,  either  theoretically  or  experimentally  up  to  the  present 
time,  w  (In  connection  with  this,  see  Supplement  7.  ) 

29.  Some  Results  of  Experimental  Investigations 


The  properties  of  medium/long  waves  differ  in  many  aspects 
from  the  properties  of  waves  10-25  km  long  and  longer.  For  ex¬ 
ample,  at  wavelengths  from  2-3  km  long  and  somewhat  longer,  the 
direct  wave,  almost  in  pure  form,  prevails  considerably  up  to 
distances  between  1000  to  2000  km  in  daytime.  After  sunset,  the 
intensity  of  the  field  increases  strongly  and  other  specific  peculiar¬ 
ities  of  medium  waves  appear  (see  para.  33),  At  the  same  time, 
the  expected  normal  diurnal  variation  of  average  values  of  the  field 
strength  of  waves  between  15-20  km  in  length  is  frequently  absent. 
Therefore,  in  the  following  discussion  of  fundamental  peculiarities 
of  long  radio  waves,  their  differences  at  the  various  frequencies  in 
this  range  should  be  kept  in  mind. 


1.  Dependence  of  Tension  of  the  Field  on  Distance 

Even  at  short  distances  from  the  radiator,  both  in  the  daytime 
and  at  night,  the  interfering  character  of  the  field,  depending  upon 
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distance,  has  an  effect  on  long  waves. 

Figure  233  shows  the  results  of  measurements  at  higher  fre¬ 
quencies;  here,  at  a  frequency  of  85  kc/s  (^-  ®  3  500  m).  The  quasi¬ 
periodical  oscillations  of  the  field  become  observable  at  distances 
from  the  radiator  greater  than  100  wavelengths,  and  increase  with 
distance.  These  variations  of  waves  can  be  interpreted  as  a  result 
of  the  superposition  of  the  direct  and  once -reflected  wave  from  the 
ionosphere.  Indeed,  the  coefficient  of  reflection  increases  with  the 
angle  of  incidence  of  a  wave  on  the  ionosphere,  which  leads  to 
strengthening  of  the  influence  of  the  reflected  wave.  Moreover, 
analysis  of  the  results  of  such  measurements  shows  that  the 
oscillations  of  the  field  frequently  occur  close  to  the  curve  of  the 
ground  wave.  At  low  and  super-low  frequencies,  the  oscillations  of 
the  field  are  greater  and  become  already  observable  at  a  distance  of 
several  wavelengths  from  the  radiator.  Thus  in  measurements  made 
at  a  frequency  of  16  kc/s  *  18,  700  km)  (Fig.  234),  the  amplitude  of 
the  field  changed  by  2  5  percent  even  at  a  5-6  k  distance  from  the 
radiator;  in  the  preceding  figure,  oscillations  of  approximately  10 
percent  began  to  be  noticeable  only  at  a  distance  where  r  »  70  k  , 
and  oscillations  of  25  percent  only  at  the  point  where  r  *  110  k.  At 
distances  between  100-200  km  and  more,  at  a  frequency  of  16  kc/s, 
it  was  impossible  to  interpret  the  results  of  the  measurements  ex-~ 
actly  by  introducing  the  direct  and  reflected  waves  into  the  analysis. 
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Fig.  233  -  Product  of  field  strength  with  distance  from  radiator. 

(Measuren'ent  on  frequency  85  kc/s  (k  ~  3  500  m) 

Translator' s  Note:  Henceforward,  many  figures  and  diagrams  in  the 
text  are  supplied  with  Russian  inscriptions  along  the  curves  and  mar¬ 
gins.  Not  being  able  to  reproduce  the  figures  in  full,  this  translator 
puts  them  in  English  in  the  same  order  as  they  appear  on  the  original, 
from  left  to  right.  The  explanatory  legend  will  be  translated  after 
the  number  of  the  figure. 
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Fig.  234  -  Dependence  of  (E  •  r)  with  distance.  Measurement  on 
frequency  16  kc/s  ~  18700  m) 


With  increased  distance,  the  variation  in  field  strength  grad¬ 
ually  smoothes  and  becomes  approximately  equal  at  various  fre¬ 
quencies,  the  amplitude  of  the  field  being  somewhat  larger  at  the 
lower  frequencies.  Experimental  averaged  curves  of  field  strength, 
depending  on  distance  (Figs.  235  and  236),  were  obtained  from  the 
data  of  many  measurements  [160]  conducted  at  different  periods  of 
time  over  the  sea,  up  to  a  distance  of  12,000  km.  In  Fig.  235,  the 
average  values  of  the  field  strength  obtained  for  day  and  mixed 
routes,  at  a  frequency  of  18  kc/s,  are  compared  with  the  formula  of 
ideal  radio  transmission,  300  \/W/r,  and  with  the  curve,  E(r),  com¬ 
puted  for  the  spherical  earth.  It  is  evident  that  as  distance  increases, 
the  experimental  values  of  E  exceed  the  values  of  E  in  the  direct 
wave  by  two  to  three  magnitudes  or  more. 

The  same  discrepancy  is  observed  for  routes  situated  com¬ 
pletely  in  the  region  of  darkness;  the  corresponding  results  of 
measurements  (Fig.  236)  are  averaged  over  the  range,  12.  5-30  kc/s. 
From  Fig.  235,  it  is  evident  that  the  amnlitude  of  the  field  may 
change  several  times  in  various  cases,  and  that  the  changes  do  not 
depend  on  frequency. 


160 

Kitchen,  F.A.  ,  D.  G.  Pressey,  andK.  W.  Tremellen,  A  review  of 
present  knowledge  of  the  ionospheric  propagation  of  very  low,  low- 
and  medium  waves,  J.  Inst.  Elec.  Eng.,  Pt.  Ill,  100  (1953). 
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showing  dependence  with  distance 
on  frequency  18  kc/s  (solid  line). 
By  dotted  line  is  shown  the  course 
of  theoretical  curves.  (W-  1  kw) 


and  range  (diapasone)  12.  5-30 
kc/s  (solid  line). 


Until  recently,  various  types  of  empirical  formulae,  obtained 
as  the  result  of  processing  a  series  of  numerous  measurements, 
were  still  used  frequently  in  the  computation  of  tension  of  the  field. 
One  such  is  the  formula  of  Austin-Cohen 


0. 0015 


_  300  nTW  /  e  ■  ^/\ 

~  r  sin  6  * 


or  Espenschied,  Anderson  and  Bailey 

0.  005 


300  nTW  I  0  \  1.25 

sin  5 


1  < 

E  =  -  /  — : — TT  e 

r  SI 


(29.1) 


(29.2) 
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There  are  supposed  to  be  present  in  the  exponents  of  the 
formulae  of  type  (29.1),  multipliers  (-  0.0014  or 

(0.0045  r/X.^*  ),  which  were  selected  each  time  to  satisfy  the  series 
of  experiments  under  scrutiny  in  the  paper. 

Other  types  of  empirical  formulae  were  also  used  whose 
derivation  is  based  on  qualitatively  correct  perception  of  the  form¬ 
ation  of  cylindrical  waves,  directed  by  the  earth  and  ionosphere,  in 
the  earth -ionosphere  waveguide.  Starting  from  this  point,  and  taking 
into  account  the  sphericity  of  the  earth,  the  flow  of  energy  from  the 
sources  crosses  the  lateral  surface 


S  =  TT  sin  6  {ZK  z  +  z®) 
o  o 


(29.3) 


237),  differing  from  the  lateral  surface,  S  =  2tt  r  z,  in  the  case 
of  a  flat  earth.  Not  taking  into  account  losses  in  the  walls  of  the 
waveguide,  we  have 


(29.4) 


and  for  the  vertical  dipole  (its  intensification  is  equal  3/2) 

tr  300  n/1^  i  r~  , 

^  =  TTTT  ThTB  (29.  5) 

where  r  and  z  are  expressed  in  km;  W  in  kw;  and  <<  1. 


Fig.  237  -  Illustrating  the  deduction  of  the 
formula  taking  into  considera¬ 
tion  the  sphericity  of  the  earth. 
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In  order  to  account  for  the  losses  in  formula  (^9.  5),  various  authors 
have  used  exponential  multipliers  of  the  type,  exp(-  a  r/^  m)  and 
have  selected  values  of  a  and  m  which  would  best  satisfy  the  results 
of  the  experiments. 

Naturally,  at  present,  there  is  no  need  to  use  formulae  of  this 
type  since  strict  theoretical  calculations  result  in  formulae  which 
permit  one  to  compute  not  only  the  field  strength  of  long  and  super - 
long  waves,  but  also  their  phase  structure  for  a  wide  range  of  fre¬ 
quencies,  and  their  dependence  on  the  conductivity  of  the  ionosphere 
and  the  earth  (see  para.  30). 

2.  Various  Properties  of  Long  Waves 

As  we  have  seen,  the  ionosphere  contributes  to  a  considerable 
increase  of  intensity  of  long  radio  waves  over  great  distances,  and 
to  the  interference  characteristics  of  the  field.  Let  us  investigate 
the  other  properties  of  long  radio  waves,  also  conditioned  by  the  in¬ 
fluence  of  the  ionosphere,  mainly  by  the  changeability  of  its  state. 

At  distances  near  the  rad'ator,  where  the  influence  of  the  direct 
wave  is  still  felt,  even  at  the  lowest  frequencies  of  this  range,  it  is 
possible  to  consider  the  field  as  a  superposition  of  direct  and  reflect¬ 
ed  waves  from  the  ionosphere.  Thus,  by  compensating  for  the  direct 
wave  in  the  experiments,  the  ionuspherically  reflected  wave  may  be 
separated.  This  has  permitted  the  determination  in  various  experi¬ 
ments  of  the  state  of  polarization  of  long  waves,  their  apparent 
height,  and  their  coefficient  of  reflection  at  vertical  and  oblique  in¬ 
cidence  [l6l] . 

Summarizing,  it  was  shown  by  various  methods,  (the  experiments 
were  conducted  within  the  range  ^  *  2-18,  7  km),  that  up  to  distances 
between  100-300  km  at  middle  latitudes,  the  reflected  wave  is  almost 
circularly  polarized.  At  greater  distances,  the  state  of  polarization 
is  changed.  For  instance,  at  ^  «  18  (7  km),  and  r  *  400  km,  purely 
linear  polarized  waves  are  already  observed.  In  this  case,  the  plane 


161 

Bracewell,  R.  C.  ,  K.  G.  Budden,  J.  A.  Ratcliffe,  T.W.  Straker  and 
K.  Weekes,  The  ionospheric  propagation  of  low-  and  very-low  fre¬ 
quency  radio  waves  over  distances  less  than  1000  km,  J.  Inst.  Elec. 
Eng.  ,  Pt.  m,  98,  221  (1951). 
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of  polarization  has  turned  approximately  45°  relative  to  the  plane  of 
incidence  of  the  wave. 

The  coefficient  of  reflection,  p  ,  of  these  waves  also  changes 
considerably  with  distance.  Figure  238  shows  the  experimental  de¬ 
pendence  of  the  coefficient  of  reflection  on  frequency,  obtained  under 
various  conditions,  almost  at  vertical  incidence;  Figs.  239  and  240 
give  the  idea  of  its  diurnal  and  seasonal  changeability  at  various  fre¬ 
quencies.  With  each  increase  in  distance,  the  coefficient  of  re¬ 
flection  increases  strongly.  For  instance,  16  kc/s  at  r  «  300  km  in 
the  described  experiments  gave  a  measured  value  of  p  ®  0.15,  and 
at  the  distances,  r  *  500-800  km,  a  value  of  p  ®  0.  33.  In  this  ex¬ 
periment,  the  coefficient  of  reflection  was  normally  observed  to 
increase  suddenly  at  a  distance  of  about  400  km.  In  summer,  at 
frequencies  of  85  and  12?  km  kc/s, the  value  of  p  ®  0.001  at  vertical 
incidence:  when  r  «  400  km,  p  »  0.01,  when  r  *  700  km,  p  *  0.04, 
and  when  r  *  900  km,  p  ~  0.  055.  In  lin^  with  this,  a  rapid  decrease 
of  p  is  observed,  over  long  distances,  one  hour  before  sunrise,  and 
a  rapid  growth  of  p  in  the  period  of  sunset.  The  apparent  altitude 
of  reflection  with  this  range  of  frequencies  changes  within  limits  be¬ 
tween  67-85  km.  At  the  lower  frequencies,  altitude  increases  from 
10-16  km  when  passing  from  day  to  night.  These  data  result  from 
observations  at  close  and  vertical  incidence  of  the  wave.  At  oblique 
incidence,  at  16  kc/s,  the  same  diurnal  variations  of  altitude  are  ob¬ 
served.  However,  at  frequencies  of  100  kc/s  and  more,  the  diurnal 
changes  in  apparent  altitude,  instead  of  being  on  the  order  of  7-10 
km,  are  about  18  km.  The  character  of  modification  of  all  the  given 
values  does  not  lend  itself  to  a  simple  explanation.  It  suggests  that 
the  lower  part  of  the  ionosphere  (its  base),  while  playing  a  big  part 
in  the  propagation  of  long  radio  waves,  has  a  series  of  peculiarities 
which  have  been  little  studied  up  to  now.  This,  in  great  part,  is  be¬ 
cause  up  to  now  there  were  no  direct  measurements  made  of  dis¬ 
tribution  of  charges  and  number  of  collisions  at  these  altitudes. 

Now,  let  us  pass  to  the  study  of  other  properties  of  long  radio 
waves.  The  state  of  the  field  of  long  waves  is  subject  to  chaotic 
variations  (so-called  fadings  or  black-outs),  as  well  as  to  regular 
changes  (diurnal,  seasonal,  and  eleven-year  cycle).  Definite  effects 
are  also  observed,  conditioned  by  variations  of  the  magnetic  field  of 
the  earth,  sudden  atmospheric  disturbances,  etc.  All  of  these  effects 
are  studied  most  completely  as  a  part  of  measurements  of  the  tension 
of  the  field. 
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Oscillations  in  amplitude  which  are  relatively  slow  a,nd  shallow 
gradually  increase  with  the  shortening  of  wavelength  and  become 
'ery  noticeable,  especially  in  the  range  of  medium  waves.  Similar 
oscillations  are  observed  in  the  recording  of  phase  and  direction  of 
arrival  of  radio  bearing  waves  (see,  for  instance,  Fig.  241). 


Local  time,  hour 


Fig.  241  -  Oscillations  of  radio  bearing  on  wave  4700  m  on  the 
distance  837  km  from  radiator. 


Let  us  note  here,  that  according  to  a  series  of  data  for  distances 
between  5000-6000  km,  the  phase  of  the  received  wave  is  almost 
constant.  This,  seemingly,  points  to  the  absence  at  this  range  of 
any  influence  by  many  modes. 

In  many  cases,  however,  the  diurnal  variation  in  amplitude  is 
subject  to  many  changes.  Figure  242  shows,  for  comparison,  the 
curves  of  diurnal  variation  of  average  monthly  values  of  field  tension, 
E,  tor  three  stations  20  kwt  strong  (k  =  5270,  11,680,  and  17,  500  m) 
over  a  distance  of  5482  km. 

From  this  figure  (drawing),  it  is  evident  that  with  the  increase 
in  wavelength,  the  limits  of  change  of  E  are  decreased,  and  that  at 
all  wavelengths,  the  value  of  E  increases  at  nighttime.  However,  it 
is  knovm  that  over  some  traverses  a  decrease  of  E  at  nighttime  is 
regularly  observed,  not  an  increase.  The  characteristic  peculiarity 
of  the  diurnal  variation  in  the  field  amplitude  of  long  waves  is,  above 
all,  the  considerable  di^crease  in  E  during  the  period  of  sunrise. 

This  phenomenon  is  established  over  a  wide  range  of  wavelengths, 
beginning  as  low  as  those  at  several  hundred  kilometers. 
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Fig.  242  -  Average  monthly  diurnal  curves  of  field  strength  on 
three  waves  at  distance  5482  km 


In  a  number  of  cases,  at  the  moment  of  sunrise,  the  amplitude 
of  the  field  has  a  deep  minimum  with  following  maximum.  After  this, 
a  relatively  constant  value  of  £  is  established,  which,  in  its  depend¬ 
ence  on  distance,  or  on  wavelength,  may  be  large  or  smaller  than 
the  night  value  of  field  amplitude.  A  similar  phenomenon  is  described 
in  the  literature  as  "twilight  effect",  because  it  is  frequently  observed 
during  sunset.  As  a  rule,  it  is  less  well-marked  at  this  time.  Twi¬ 
light  effect  is  absent  altogether,  both  in  the  morning  and  evening,  on 
many  strip-chart  recordings  (Fig.  243).  The  seasonal  course  of  field 
amplitude  is  also  subject  to  changes.  Amplitude  of  field  is  greater 
in  daytime  over  many  paths,  and  greater  in  summer  than  in  winter; 
at  nighttime,  it  is  approximately  the  same  all  year  round. 
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Fig.  243  -  Diurnal  course  of  field  amplitude 
a  -  on  frequency  16  kc/s  on  distance  2500  km; 
b  -  on  distance  17.  8  kc/s  on  distance  7860  km; 
c  -  on  frequency  26  kc/s  on  distance  11200  km. 


There  are  some  indications  that  the  amplitude  of  field  is  greater 
over  paths  situated  at  the  same  latitudes  than  over  paths  along  the 
same  meridians. 

All  these  effects  are  explained  by  the  complicated  structure  of 
the  field  of  long  radio  waves,  depending  upon  distance  from  the 
radiator,  degree  of  illumination  of  route,  height  of  the  base  of  the 
ionosphere,  character  of  transition  of  altitude  into  region  of  semi¬ 
darkness,  and  similar  phenomena. 
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The  periods  of  sudden  disturbances  in  the  ionosphere  are 
characterized  by  the  fact  that  at  this  time  there  occurs  a  pronounced 
weakening  of  the  field  of  short  waves;  the  field  strength  of  the  long 
waves,  on  the  contrary, often  increases  exactly  at  this  time.  This 
can  be  seen  in  Fig.  244,  which  shows  an  increase  of  the  field  ampli¬ 
tude  at  a  wavelength  of  11,000  meters  during  a  period  of  sudden  dis¬ 
turbance  which  began  at  11  h  57  m.  At  the  time  of  sudden  disturbances, 
as  we  have  seen,  the  structure  of  the  E-  and  R-regions  almost  does 
not  change;  only  their  electronic  and  ionic  concentrations  increase, 
with  a  consequent  increase  in  the  gradient  of  ionization  of  the 
D-region.  Therefore,  at  this  time,  the  absorption  of  waves  passing 
through  the  D-layer  increases,  and  the  losses  in  intensity  of  long 
waves  directed  by  this  layer  may  decrease  because  they  leak  into 
the  D-layer  less.  With  an  increase  of  frequency  in  tliis  range,  a 
decrease  ,  and  not  an  increase,  in  the  amplitude  of  the  field  may 
well  prevail  at  the  time  of  sudden  disturbances. 


Fig.  244  -  The  growth  of  the  amplitude  of  field  on  wave  11100  m 
in  the  period  of  sudden  disturbance  which  began  at  11  h  57  m. 


Numerous  observations  have  shown  that  with  the  increase  in 
activity  of  the  sun  and  magnetic  storms,  the  field  strength  of  the 
long  waves  (Fig.  245)  increases;  the  reverse  is  observed  for  medium 
waves.  Such  dependence  is  also  apparently  explained  by  the  increase 
in  the  gradient  of  electronic  concentration  in  the  base  of  the  ionosphere 
during  these  periods. 
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In  general,  computations  of  the  field  of  long  radio  waves  de¬ 
mand  a  solution  of  the  waveguiding  problem  for  a  spherical  earth, 
surrounded  by  a  nonhomogeneous  ionosphere,  the  properties  of 
which  are  modified  substantially  with  frequency.  In  such  a  situation, 
it  is  difficult  to  solve  this  problem  mathematically.  Up  to  now,  no 
one  had  succeeded  in  reducing  the  problem  to  simple  form  which 
would  permit  a  discussion  of  the  experimental  data  on  this  basis. 

At  the  same  time,  a  great  number  of  papers  were  devoted  to  this 
problem  [162  -  164],  In  the  latest  series  of  investigations 

162 

Watson,  G.  N.  ,  The  transmission  of  electric  waves  around  the 
earth,  Proc.  Roy.  Soc.  ,  _95,  546  (1919). 

163 

Schumann,  W.  ,  Uber  die  Oberfelder  bei  der  Ausbreitung  langer 
elektrische  Wellen  in  System  Erde-Luft-Ionosphare,  Zs,  f.  angew. 
Phys.,  4,  474  (1952);  6,  35,225,  267,  and  346  (1954). 

164 

Krasnushkin,  P.  E.  ,  N.  A.  Yablochkin,  The  theory  of  propagation  of 
very  long  waves  applied  to  the  problem  of  long  distance  communi¬ 
cation  and  navigation,  Proc.  of  the  State  Union,  Scientific  Research 
Institute  (1955). 
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[163,  165],  such  a  spherical  problem  was  formulated  for  a  homo¬ 
geneous  ionosphere.  In  the  process  of  computation,  however,  the 
authors  preserve  the  influence  of  the  sphericity  of  the  earth  only  by 
d/sin  6  .  Thus,  as  special  cases,  they  obtain  previously  investigated 
solutions  for  the  plane  problem  with  a  uniform  upper  boundary  [166, 
167].  In  these  papers,  allowing  for  sphericity  is  reduced  to  taking 
into  account  a  ’’geometric”  factor,  also  obtainable  from  simple  con¬ 
siderations  (see  deduction  (29.  5)). 

Consequently,  the  spherical  problem  for  a  homogeneous 
ionosphere  is  solved  in  one  of  the  most  recent  papers  [169]  for 
sufficiently  high  conductivity  of  earth  when  the  approximate  boundary 
conditions  are  suitable.  The  equation  obtained  in  this  paper,  de¬ 
fining  the  wave  numbers  of  wave  spectrum,  is  rather  complicated. 
The  author  does  not  solve  it.  However,  using  this  equation  approxi- 
nnately,  it  is  possible  to  correct  the  solution  for  the  flat  earth  and, 
consequently,  take  into  account  the  influence  of  its  sphericity.  (See 
Supplement  6  in  which  results  of  this  work,  and  also  of  paper  [242] , 
published  after  this  present  work  went  to  print,  are  given  briefly. 
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Schumann,  W.  ,  Uber  die  Oberfelder  bei  der  Ausbreitung  langer 
elektrische  Wellen  in  System  Erde-Luft-Ionosphare. 

165 

Wait,  J.R.  ,  The  mode  theory  of  VLF  ionospheric  propagation  for 
finite  ground  conductivity,  Proc.IRE,  45,  760-768  (1957). 

166 

Brekhovskikh,  L.  M.  ,  Waves  in  stratified  media,  Acad.  Sci,  ,  SSR,  1957; 
Bui.  Acad.  Sci.  ,  Physics  Series,  (1949,  p.  505;  ^  pp.  515-534. 

167 

Budden,  K.  G.  .  The  propagation  of  a  radio  atmospheric,  Phil.  Mag.  , 
42,  1  (1951);  «,  1179  (1952). 

169 

Bremmer,  H.  ,  National  Bureau  of  Standards  Report  No.  5518  (1957). 
242 

Wait,  J.R.  ,  Terrestrial  propagation  of  very  low-frequency  radio 
waves,  a  theoretical  investigation,  J.  Research  NBS  64D (Radio 
Prop.),  No.  2,  153-204  (Mar. -Apr.  I960). 
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One  of  the  basic  conditions,  without  which  it  is  innpossible  to 
solve  the  waveguiding  problem  for  long  radio  waves  correctly,  is 
the  need  to  take  into  account  the  nonhomogeneity  and  dispersion  of 
the  ionosphere,  and  the  dependence  of  its  reflective  properties  on 
frequency  which  is  modified  considerably  in  the  frequency  range 
under  consideration.  It  is  very  difficult  to  solve  such  a  problem 
mathematically  in  actual  form;  all  the  more  since  the  analytical 
formulae  describing  the  modification  of  corresponding  parameters 
of  the  ionosphere  are  not  known.  Therefore,  one  is  forced  to  use 
methods  of  computation  based  on  the  physical  picture  of  phenomena, 
and  to  describe  the  real  properties  of  the  ionosphere  with  the  aid  of 
effective  parameters  and  equivalent  parameters  of  a  homogeneous 
ionosphere.  Such  a  method  of  computation  for  a  flat  earth  is  laid 
down  in  paper  [168].  As  analysis  has  shown,  it  can  be  reconciled 
very  well  with  v;xperimental  data  over  a  wide  range  of  frequencies. 

The  method  of  computation  is  given  briefly  below,  together  with  re¬ 
sults  of  some  comparisons  between  numerical  computations  and 
experimental  data. 

1.  Method  of  Computation 

The  method  for  nonhomogeneous  plane  waves,  used  previously 
for  computing  the  field  of  the  dipole  over  the  flat  earth  [170],  is 
added  to  the  basic  steps  for  solution  of  the  flat  waveguiding  problem. 
This  method  has  an  obvious  physical  basis.  Let  us  describe  it 
briefly,  since  this  permits  us  to  understand  how  the  solution  for  the 
homogeneous,  sharply  bounded  ionosphere  is  transferred  directly  to 
the  case  of  a  nonhomogeneous  (stratified)  ionosphere. 

As  is  known,  the  field  of  an  elementary  dipole,  proportional  in 
free  space  to  exp(ikR)^,  can  be  resolved  into  a  continuum  of  non¬ 
homogeneous  plane  waves.  Therefore,  formally,  the  field  of  a  dipole, 
situated  over  or  near  the  boundary  of  separation,  may  be  described 
as  a  result  of  superposition  of  a  direct  wave  and  of  a  continuum  of  non- 
homogeneo  s  plane  waves  reflected  from  this  boundary.  The  co¬ 
efficient  of  reflection  from  the  separation  boundary  is  generalized  in 


168 

Al'pert,  Ya.  L.  ,  On  the  propagation  of  low  frequency  electromagnetic 
waves  over  the  earth' s  surface,  Bui.  Acad.  Sci.  ,  USSR,  Phy.  Series, 
1955;  Progress  in  Physics,  ^(1956),  p.  369. 

170 

Weyl,  H.  ,  Ausbreitung  elektromagnetischer  Wellen  uber  einen 
ebenen  Leiter,  Ann.  d.  Phys.  ,  60,  481  (1919). 
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this  case  in  such  a  way  that  the  sines  and  cosines  of  angle  of  in¬ 
cidence  are  substituted  for  by  their  complex  values.  These,  in 
turn,  serve  as  arguments  of  plane  nonhomogeneous  waves.  The 
corresponding  consideration,  applied  to  a  plane  waveguide,  leads  to 
the  appearance  of  four  groups  of  plane  nonhomogeneous  waves  [166  - 
168]  whose  superposition  gives  the  following  values  of  field  when  the 
source  and  the  point  of  observation  are  located  on  the  boundary  of 
the  separation:  Vertical  and  horizontal  components  of  electric  and 
magnetic  fields 
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where  E  and  H  are  expressed  in  mv/m;  where  h  is  the  effective 
height  (of  the  transmitting  antenna),  z  is  the  height  of  the  ionosphere; 
^  is  the  wavelength;  r  is  the  distance  in  km;  and  I  is  the  current  in 
amperes.  In  formula  (30.1),  H'  ^  and  are  Hankel' s  functions 
with  the  function,  P(C  )  given  by 
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Budden,  K.  G.  ,  The  propagation  of  a  radio  atmospheric. 
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waves  over  the  earth' s  surface. 
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Pj(C)  and  Pg  (C)  are  formulae  of  coefficient  of  reflection  from  the 
surface  of  the  earth  and  the  ionosphere.^  The  complex  numbers, 

(S^  and  C^),  expressed  as 

S  =nA-C2=S  +S  ,  (30.3) 

n  n  nj  ng 

are  sine  and  cosine,  and  define  the  wave  numbers,  k  S  of  the 
spectrum  of  cylindrical  waves,  excited  in  the  earth -ionosphere 
waveguide.  They  are  computed  from  the  equation  of  poles,  obtained 
from  the  solution  of  this  problem  and  having  the  form 

2h  k  C  +  i  In  [p,(C  )  Pg(C  )]  =2n  TT  .  (30.4) 

on  ^  n  *  n 


We  now  use  the  asymptotic  forms  of  the  Hankel  functions,  which 
are  well  satisfied,  if 

Ik  S  r|>4~5,  2n<  arg  (k  S  r)  <  rr  ,  (30.  5) 

on  on 

then  introducing  the  multiplier  300  <s/W  instead  of  120  nl  h,  and  also 
the  geometric  factor,  which  takes  into  consideration  the  sphericity 
of  the  earth,  we  obtain 
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The  factor  l/2  in  eqn.  (30.  2)  should  be  l/4.  It  should  be  noted  that 
this  factor  was  only  given  in  the  author' s  previous  work  [168]  for 
pj  =  1.  Its  generalization  appears  in  reference  [242],  (JRW). 
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Thus,  it  is  evident  that  the  field  is  a  superposition  of  the  spectrum 
of  cylindrical  waves,  decreasing  in  a  proportion  of  the  form 
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In  formula  (30.  6)  the  function  standing  under  the  sum  is  an 
analogue  of  the  attenuation  function  f  (p  ),  but  in  contrast  to  f(p ),  its 
modulus  does  not  decrease  monotonically  with  distance,  but  oscillates. 
Therefore,  the  formula 
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(30.7) 


may  be  called  an  interference  multiplier  of  the  field  of  radio  waves 
in  the  waveguide.^ 

The  formulae  given  above  show  that,  for  the  computation  of  field 
of  long  waves,  it  is  necessary  to  determine  the  wave  numbers  k  S  . 
The.  e  are  obtained  by  solving  a  rather  complicated  transcendental 
equation  of  the  poles  (30.4),  Here  Pi(C)  is  the  Fresnel  coefficient 
from  the  earth.  As  for  the  form  of  function,  pg  (C),  it  is  unknown, 
inasmuch  as  the  coefficient  of  reflection  from  the  nonhomogeneous 
ionosphere,  generally  speaking,  modifies  the  various  boundary  con¬ 
ditions.  Having  this  difficulty  in  view,  it  is  possible  to  solve  this 
problem  by  the  following:  The  equation  (30.4)  is  first  computed  for 
a  sufficient  range  of  electric  parameters  of  the  ionosphere  and 
various  frequencies,  selecting  as  p,,(C)  the  Fresnel  coefficient  of 
reflection  from  the  homogeneous  medium.  This  permits  us  to  obtain 
the  dependence  of  poles,  on  these  parameters.  Then  we  compute 

for  various  frequencies,  and  for  a  given  real  model  of  ionosphere, 
the  true  values  of  modulus  and  of  the  argument  of  reflection  coefficient. 
Next,  after  comparison  with  the  corresponding  values  of  Fresnel' s 
coefficient  of  reflection,  we  compute  the  dependence  of  the  medium 

r 

Apparently  in  writing  eqn.  (3C.  7),  P(C  )  has  been  set  e:  al  to 

one.  (JRW)  " 
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on  frequency  and  electric  parameters.  Thus,  the  effective  values  of 
corresponding  parameters  give,  with  the  help  of  Fresnel' s  formula, 
the  same  values  of  intensity  and  phase  for  the  reflected  wave.  The 
effective  value  of  the  coefficient  of  refraction  found  by  this  method  is 
used  later  on  for  determination  of  and  This  method  was  follow¬ 

ed  in  paper  [168] . 


Below  are  the  corresponding  results  of  numerical  computations 
for  vertically  polarized  waves.  The  generalized  Fresnel  formula  has 
the  form; 
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neglecting  the  influence  of  the  magnetic  field  of  the  earth, 
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3 
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4tt  Ne^ 

m  V®  y 


4ne®  N_ 
m  o)  V  ‘ 


Then,  if 
(30.9) 


Since  the  first  term  in  (30.9)  is  near  to  one,  then  for  the  lower  part 
of  the  D-region,  which  plays  the  basic  part  in  transmission  of  long 
radio  waves  in  daytime,  it  is  suitable  for  this  time  of  day  to  select 
relation  (N/v)g££  as  the  effective  parameter  of  medium  [168],  defined 
by  adjusting  Fresnel' s  coefficient  to  the  true  coefficient  of  reflection. 
Therefore,  if  the  theoretical  dependences  of  pg(C)  and  on  (N/v)are 
established,  then  the  first  of  them  determines  (N/vg££),  and  the 
second,  the  corresponding  value  of  to  it. 

2.  Dependence  of  Field  Strength  Distance  and  Frequency 


We  have  seen  that  the  structure  of  the  field  of  long  waves  can  be 
described  by  interference  multiplier  (30,  7).  The  analysis  of  the 
modulus  (30.  7),  shows  that  the  amplitude  of  field,  B  (  u»,  r),  changes 
in  an  irregular  way,  depending  upon  distance  from  its  source,  since 
the  spatial  periods,  or  the  lengths  of  waves,  of  various  com¬ 

ponents  of  the  spectrum  are  comparable.  On  moving  away  from  the 
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source,  the  factor  B(c)),  r)  gradually  becomes  qua  si -periodic  since 
the  majority  of  spectrum  waves  (modes)  are  damped  with  the  ex¬ 
ception  of  the  first  two  or  three  (Fig.  246).  The  amplitude  of  field 
changes  in  a  number  of  details  with  relatively  small  changes  in 
ionospheric  parameters.  This  is  especially  true  near  the  radiator 
(Fig.  246  a,  b,  c. ).  The  same  phenomenon  is  observed  with  small 
changes  of  frequency  in  constant  ionospheric  parameters  (Fig.  247). 


Fig.  246  -  Dependence  of  modulus  of  interference  multiplier  B(w,  r) 
from  distance  for  various  values  of  N/v  on  frequency  of  15  kc/s. 


Fig.  247  -  Dependence  B(u,  r)'s/r  on  distance  on  two  near  frequencies. 
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The  results  of  computations,  conducted  for  real  modes  of  the 
ionosphere,  showed  that,  in  daytime,  beginning  with  30  kc/s  and  be¬ 
low,  and  over  distances  exceeding  1000  km,  the  basic  part  is  played 
only  by  the  senior  component  (dominant  mode)  of  the  wave  spectrum 
(so-called  zero  mode).  At  nighttime,  the  number  of  important  modes 
increases.  For  example,  at  f  *  15  kc/s,  it  is  necessary  to  take  into 
account  more  than  5-7  modes;  at  a  frequency  of  30  kc/s,  more  than 
10-12;  there  are  even  3  or  4  modes  at  f  *  5  kc/s. 

The  necessity  of  taking  into  account  the  large  number  of  terms 
of  the  sum  which  determines  the  interference  multiplier  makes  com¬ 
putation  of  the  field  for  nighttime  more  difficult.  Another  difficulty 
lies  in  the  fact  that  it  is  more  complicated  to  compute  the  true  co¬ 
efficient  of  reflection  from  the  ionosphere  for  nighttime.  This  is 
conditioned  by  the  fact  that  in  the  daytime  the  coefficient  of  refraction 
has  a  simpler  form  and  depends  only  on  one  parameter,  N/v,  so  that 
it  can  be  described  in  the  range  of  low  frequencies,  with  which  we 
have  been  dealing,  with  the  aid  of  formula  (30.9).  The  course  of 
N/v  under  various  conditions  is  close  to  the  curve  shown  in  Fig.  248. 
Therefore,  as  long  as  the  frequency  is  not  excessively  small  (so  that 
the  length  of  the  wave  is  on  the  order  of,  or  smaller  than,  200-250  km, 
--  the  thickness  of  the  lower  half  of  the  ionosphere),  and  the  leakage 
through  the  ionosphere  of  waves  at  low  frequency  is  not  large,  it  is 
possible  to  approximate  the  ionosphere  by  a  transition  layer.  This 
assumes  (N/v)  *  a  constant  greater  than  89-90  km.  This  reduces 
computation  considerably. 


Z.Kf^ 


Fig.  248  -  Dependence  of  n/v  on  the  altitude  in  the  lower  part  of  the 
ionosphere  in  daytime.  The  curve  is  drawn  which  closely  approximates 
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Curves  of  the  field  of  strength,  as  a  function  of  frequency  for 
various  distances  (Figs.  249-252),  were  computed  with  the  aid  of 
such  a  model  of  the  ionosphere  by  the  method  described  above.  In 
computations,  it  was  assumed  that  the  coefficient  of  reflection  from 
the  earth  is  equal  to  one,  i.  e.  ,  o-  =  oo.  The  values  of  and  P(C„), 
obtained  for  the  given  model  of  the  ionosphere,  are  given  in  Table  37, 
as  are  also  values  of  (N/v)g££,  which  characterize  the  effective  con¬ 
ductivity  of  the  lower  region  of  the  ionosphere  at  various  frequencies. 
In  computation,  a  number  of  modes  were  used  whose  corresponding 
values,  and  P(C^),  are  given  in  Table  37. 


Frequency,  kc/s 

Fig.  249  -  Curves  of  field  strength  at  daytime  in  the  range  500-30,000 
kc/s  on  distances  100,  500  and  1000  km. 
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Table  37 

Values  Sj^,  P(C  )  and  (N/v)gff  of  various  modes  for  the 
daytime  modeTof  the  ionosphere. 
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Fig.  250  -  Curves  of  the  field  strength  ^51  -  Curves  of  the  field  strength  in 

in  the  range  500-30,  000  c/s  on  distances  range  500-30,000  c/s  on  distance 

1000,  2000  and  3000  km.  5qqq 


Distance,  km 

Fig.  252  -  Curves  of  the  field  strength  in  the  range  50-30,000 
c/s  on  distance  of  10,  000  km. 


Distance,  km 

Fig.  253  -  Results  of  measurements  (a)  field  strength  on  16 
kc/s,  measured  in  an  aircraft  over  England. 


Distance,  km 

Fig.  254  -  Results  of  measurements  (a)  and  computations  (b) 
of  the  field  strength  on  the  frequency  16  kc/s, 
measured  in  an  aircraft  on  the  route  England-Cairo 
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Examining  Figs.  249-252,  which  show  only  average  curve  of 
field  strength  (the  separate  oscillations  are  smoothed  out),  indicates 
that  the  curve  of  field  strength,  as  a  function  of  frequency,  is 
characterized  by  the  following  peculiarities:  For  a  constant  intensity 
of  the  source  (as  a  function  of  frequency),  over  the  range  7-14  kc/s, 
the  amplitude  has  a  smooth  maximum  for  all  distances.  After  this,  it 
slowly  decreases.  With  increase  in  distance,  the  maxima  are  dis¬ 
placed  towards  greater  frequencies.  A  sharp  minimum  of  amplitude 
of  the  field  occurs  at  frequencies  between  2-3  kc/s.  Its  depth  in¬ 
creases  strongly  with  distance,  so  that 

E  E 

»  10  at  r  «  1000  km,  and  »  2  •  10^ at  r  «  3000  km. 

min  min 


Besides,  there  is  a  second  maximum  amplitude  in  the  ultralong  fre¬ 
quency  legion,  as  shown  in  the  following  paragraph. 

Comparing  the  results  of  various  measurements  with  compu¬ 
tations  gives  rather  good  agreement,  if  the  given  model  of  the 
ionosphere  is  used,  even  without  detailed  precise  definition  of  the 
value  of  (N/v)  in  order  to  conform  to  conditions  of  the  experiment. 
This  is  evident  from  Figs.  253  and  254  which  show  the  results  of 
measurements  and  computations  of  the  field  tension  at  a  frequency  of 
16  kc/s.  Naturally,  for  a  closer  analysis  of  the  results  of  measure¬ 
ments,  it  is  necessary  to  know  the  dependence  of  N/v  on  altitude. 


3.  Phase  and  Velocity  of  Radio  Waves 

Analysis  of  the  argument  of  the  complex  interference  multi¬ 
plier  shows  that  phase  also  changes  in  an  irregular  way  with 
distance  near  the  source,  and  then,  asymptotically,  tends  to  the 
linear  factor  (k^  S  r),  when  all  waves  are  damped  except  the  zero 
mode.  The  corres^nding  regularities  are  evident  from  Fig.  255, 
which  shows  the  course  of  the  supplemental  phase 

(30.10) 


for  the  frequency. 


15  kc/s,  and  for  two  values  of  (N/v) 

eff 


t 

In  Fig.  255,  apparently,  he  has  plotted  -  (p  (JRW) 
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Natu  rally, 
similar  way 


the  differential  phase  velocity  also  changes  in  a 


V 


9  r 


b) 

Hr 


(30.11) 


(Fig.  256).  and  the  average  phase  velocity 


c 


u  r 


(30.12) 


(Fig.  257).  Differential  velocity  deflects  in  an  irregular  way  for 
different  distances  from  the  radiator.  It  may  be  larger  or  smaller 
than  c,  and  it  departs  from  c  by  10  to  15  percent  and  more.  On  the 
other  hand,  the  average  phase  velocity  is  always  greater  than  c  by 
3  to  4  percent  near  the  source,  and,  if  far  from  it,  only  by  a  fraction 
of  a  percent.  In  the  same  manner  as  the  amplitude  of  the  field,  the 
curve  of  the  phase  and,  consequently,  the  velocity,  changes  in  de¬ 
tail  with  small  changes  of  ionospheric  parameters  or  frequency, 
while  preserving  its  general  character  (see  Figs.  255-257). 


Using  the  sanie  model  of  the  ionosphere,  the  dependence  of 
average  phase  velocity  of  radio  waves  on  frequency  for  various 
distances  from  the  radiator  is  computed.  The  corresponding  de¬ 
pendences  are  given  in  Figs.  258-260,  Here,  considering  the  small 
difference  between  v  and  c,  values  art  derived  of  (v/c-1).  From 
these  figures,  it  is  evident  that  average  phase  velocity  is  everywhere 
larger  than  c,  and  decreases  monotonically  in  the  range  of  frequencies 
above  its  maximum,  between  2  and  3  kc/s(v-  c)  >  0.  From  the 
reasoning  in  the  following  paragraph,  where  experimental  values  (ob¬ 
tained  from  the  analysis  of  the  spectrum  of  atmospherics)  are  com¬ 
pared  to  theoretical  values  of  v/c,  it  follows,  however,  that  at 
f  ^  2  kc/s,  there  is  a  discrep>ancy.  It  is  supposed  [171]  that  this  dis- 
crep>ancy  depends  on  the  fact  that,  at  these  frequencies,  the  selected 
model  of  the  ionospheric  transition  layer  becomes  unsuitable  since  the 
wavelength  approximates  the  thickness  of  the  whole  lower  ionosphere 
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and  more.  Also,  the  influence  of  leakage  of  waves,  excited  in  the 
earth-ionosphere  waveguide,  into  outer  space  beyond  earth  limits, 
is  felt.  Such  leakage  is  already  noticeable  at  frequencies  between 
7  and  10  kc/s.  This  is  the  condition  which  obtains  in  the  case  of  the 
appearance  of  the  so-called  whistling  atmospherics,  described  in 
the  following  paragraph.  Up  to  now,  sufficiently  extended  calculations 
have  not  existed  for  the  propagation  of  electromagnetic  waves  at  low 
frequency  which  took  into  account  this  semi-transparency  of  the 
ionosphere. 


Fig.  257  -  Dependence  of  relation  of  average  phase  velocity  v  to  c 
on  distance  for  frequency  15  kc/s  for  three  different  values  of  N/v. 
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31.  Electromagnetic  Waves  of  Ultralow  Frequency 

This  range  of  electromagnetic  waves,  created  by  various  radio¬ 
arrangements,  is  limited  to  long  radio  waves  from  20  to  25  km  long 
(15-12  kc/s).  However,  a  source  exists  in  nature,  radiating  a  con¬ 
tinuous,  especially  intense  spectrum  of  waves,  beginning  with  fre¬ 
quencies  from  30  to  50  kc/s,  and  extending  down  to  several  hundred 
c/s  and  less.  This  source  is  lightning  discharges.  The  waves  ex¬ 
cited  by  them  are  received  over  long  distances.  It  is  possible  to 
record  several  such  signals  (called  atmospherics)  a  minute,  almost 
continuously  at  every  point  of  the  globe.  The  investigation  of  atmo¬ 
spherics  permits  us  to  obtain  a  series  of  data  on  the  propagation 
properties  of  the  waves  forming  a  part  of  their  spectrum.  Using 
this  very  method,  it  was  possible  to  determine  the  average  velocity 
vs.  time  of  electromagnetic  waves  in  the  range  from  3-20  kc/s  [172], 

The  analysis  of  atmospherics  has  also  permitted  us  to  show  that 
the  minimal  losses  near  the  earth  are  suffered  by  waves  between 
25-30  km,  and  to  obtain  valuable  data  on  the  properties  of  the  lower 
ionosphere  [173]  which  play  an  important  part  in  propagation  of  the 
whole  range  of  radio  waves.  Finally,  based  on  theoretical  analysis 
of  the  form  of  atmospherics,  we  can  check  to  what  degree  methods 
of  computation  of  long  waves  (currently  in  use  and  described  above), 
are  adequate  tools  for  the  treatment  of  the  physical  substance  of  the 
phenomena  over  a  wide  range  of  low  frequencies  [173].  In  line  with 
this,  experimental  investigations  of  some  varieties  of  such  signals  - 
the  so-called  whistling  atmospherics  -  taken  in  conjunction  with  the 
theory  of  their  propagation,  will  permit  us  to  estimate  the  leakage 
of  long  waves  through  the  ionosphere,  and  to  obtain  the  properties 
of  the  outer  ionosphere  and  interstellar  gas.  It  also  will  allow  us  to 
obtain  important  information  about  the  shortest  waves  in  this  range. 
Thus,  the  study  of  propagation  of  electromagnetic  waves,  excited  by 
lightning  discharges,  is  closely  linked  at  present  with  investigations 
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of  the  propagation  of  radio  waves.  Below,  in  brief,  are  given  sontie 
properties  and  results  of  investigations  of  atmospherics,  which  are 
of  interest  from  exactly  this  point  of  view. 

1.  Atmospherics 

An  atmospheric  is  a  signal,  formed  by  the  spectrum  of  radiated 
lightning,  at  a  sufficiently  large  distance  from  it.  Many  near-field 
observations  show  that  among  various  forms  of  discharges  of 
lightnings,  the  most  frequently  registered  signals  are  of  the  type 
presented  in  Fig.  261  a,b;  while  at  greater  distances,  exceeding 
200-300  km  approximately,  in  75  percent  of  the  cases,  the  received 
signals  are  of  the  type  shown  in  Fig.  261  c,  d. 


It  is  established,  that  for  various  computations,  at  frequencies 
seemingly  exceeding  1-2  kc/s,  it  is  possible  to  select  the  source  of 
atmospherics  of  more  or  less  standard  form,  described  by  a  tempo¬ 
rary  function 


E(t,0)  »  e‘ 


at  -  St 
-  e 


(31. 1) 


and,  consequently,  having  the  spectral  density 

arctg  {  -  - -t) 

_  ®  a  3  -  wP 


A(w,  0)  e 


ill'  («.  0) 


_ _ 


(31.2) 


with  the  values  of  a  *  4.  4  •  10^  and  p  »  4.  6  *10  [241] 


t 

To  be  consistent,  i  ^  (o).  0)  in  (31.  2)  should  be  prefixed  with  a 
negative  sign  (JRW) 
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Fig.  261  -  Photo-oscillograms  of  lightning  discharges  at  near 
distances  (a  and  b)  and  signals  excited  by  them  (c  and  d)  on 
large  distances  The  sinusoid  on  the  lower  part  of  c  and  d 
is  drawn  for  a  nme  scale,  (nnksec  =  micro  sec) 
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Temporary  function,  E(t,  r),  describing  the  form  of  the  atmo¬ 
spheric  itself,  is  the  Fourier's  integral 
00 

E(t,r)=,j^  f  {BU  r)  A(„,  0)}  e' ‘  e‘ “  ‘  d,.  , 

(31.3) 


-  00 


where  B(«,  r)  e^  is  the  interferential  multiplier  obtained  above 
[see  (30.  7)],  or,  alternatively  speaking,  a  function  of  the  passage  of 
an  atmospheric  through  the  waveguide. 


Naturally,  for  the  study  of  the  observed  experimental  data,  it 
is  expedient  to  approach  signals, E(t,  r)  from  two  aspects.  First,  it 
is  possible  to  set  a  general  condition  of  study  in  summary  form  over 
a  wide  range  of  frequency  for  the  propagation  function  ,B(  w,  r)  e”  ^ 

For  this,  it  is  necessary  to  perform  a  theoretical  synthesis  of  waves 
propagating  from  the  lightning  discharge  by  computing  the  integral 
(31.  3),  and  checking  to  what  extent  the  theoretical  form  of  the  signals, 
E(t,  r),  coincides  with  the  experimental  one.  Such  investigations 
should  show  the  general  deficiencies  of  the  theory. 


Secondly,  it  is  also  possible  to  solve  the  reverse  problem  on  the 
basis  of  harmonic  analysis  of  E(t,  r);  that  is,  by  determining  the  linked 
functio^  (31.  3), to  find  the  interference  multiplier  from  the  observed 
signal  ' 


B(«,  r) 

A(oj,  0) 


(31.4) 


where  oo 

..  .  iiKw.r)  C  t-/..  V  -iwt 

A(w,  r)  e  =  \  E(t,  r)  e  , 

J  dt  . 

o 

In  this  way  it  is  possible  to  study  the  properties  of  interference 
multiplier  by  the  analysis  of  atmospherics.  And,  from  the  harmonic 


should  be  e'  ^  *  (JRW) 


tt  . 

In  the  right-hand  side  of  (31.4)  e  should  be  e  and  in  the  following 
equation  e^'^  should  be  e"  *  ^  (JRW) 
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analysis  of  E(t,  r),  we  are  able  to  determine  the  phase,  4^(u,  r),  with 
precision  to  a  linear  term;  that  is,  only  the  so-called  supplementary 
phase  is  obtained.  Thus,  the  phase  characteristics,  obtained  from 
such  analysis,  determine 

(w,  r)  =  *  (o),  r)  -  —  =  r)  -  iji  (w,  0)  (31.  5) 

c 

and, consequently,  the  dependence  of  average  phase  velocity  of 
electromagnetic  waves  at  a  given  frequency  over  a  given  distance 
[172]. 

Further  analysis  of  the  results  obtained  by  the  method  described 
above  leads  to  the  discovery  of  ionospheric  parameters  which 
characterize  the  propagation  of  these  waves.  For  instance,  given 
theoretical  dependence,  the  velocity  of  radio  waves  may  be  obtained 
from  these  parameters. 

This  shows  us  that  the  proposed  course  of  investigation  has  been 
carried  out  in  full.  First,  theory  is  checked  by  comparison  with 
experiments.  Then,  physical  values  are  determined  from  the  experi¬ 
ments  which  characterize  the  propagation  of  the  waves  under  con¬ 
sideration.  Finally,  if  the  validity  of  the  theory  is  proven,  then  the 
values  obtained  by  experimental  means  are  again  compared  with  the 
theoretical  dependences  of  those  values  of  the  parameters  of  the 
medium,  and,  in  this  way,  the  latter  determined.  In  the  quoted 
papers  [168,  171-1  73], a  large  part  of  such  a  program  of  investigations 


168 

Al'pert,  Ya.  L.  ,  On  the  propagation  of  low  frequency  electromagnetic 
waves  over  the  earth' s  surface. 

171 

Al'  pert,  Ya.  L.  ,  and  S.  V,  Borodina,  An  investigation  of  long  and  very 
long  radio  wave  propagation  using  analysis  of  atmospherics. 

172 

Al'pert,  Ya.  L.  ,  and  S.  V.  Borodina,  On  the  propagation  velocity  of 
audio  frequency  electromagnetic  waves. 

173 

Al'pert,  Ya.  L.  ,  and  D.  S.  Fligel,  A  synthesis  of  atmospherics  and 
the  effective  parameters  of  the  lower  part  of  the  ionosphere  at  low 
frequencies. 
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was  completed.  The  results  obtained, which  immediately  character¬ 
ize  propagation  of  electromagnetic  waves  of  radio  range  and  ultralow 
frequency,  are  briefly  given  below. 

For  the  synthesis  of  wave  forms  of  atmospherics  in  the  range 

500-30,  000  c/s  ,  the  values  of  wave  numbers,  kg  (shown  above  in 

Table  37),  are  utilized.  At  frequencies  smaller  than  500  c/s  and 

greater  than  30,  000  c/s,  the  approximate  solution  of  poles  (30.4)  is 

useable  for  determining  k  S  . 

on 

If  5  <  f  <  500  c/s,  then  ^ 


2  \f2  zk  <v/S 
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(31.  6) 


where 


n  =  1,  2 . 


a  = 


nn 

zk 


P  = 


4TT{r  4ne®  N 


u  m  V 


Eqns.  (31.  6)  for  the  S  '  s  are  only  valid  if  z  k  ^^  »  1,  see 
reference  [242] ,  (JR'\^).  ° 


If  f  ^  5  c/s,  then^ 
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(31.  7) 


For  frequencies  where  f  >  30,000  c/s,  the  coefficient  of  reflection 
from  the  homogeneous  ionosphere  is  well  approximated  by  the  formula 


p(C)  =  e 


-(Bj  +  i  Bg)  C. 


(  31.  8) 


Therefore,  if  B^  and  Bg  are  computed  by  comparing  p  (C)  with  the 
true  coefficient  of  reflection  of  the  ionosphere,  it  is  possible  to  solve 
the  equation  of  poles  simply. 

With  the  help  of  these  formulae,  the  functions  {[B(b>,  r)  A(w,  r)]}  , 
and,  {[4*  (oj,  0)  +  4^(cj,  r)]},  (Figs.  262  and  263),  were  obtained,  and 
then,  by  numerical  integration,  (31.3),  signals,  £(t,r),  computed. 

It  should  be  kept  in  mind  when  evaluating  the  curves  in  Figs.  262 
and  263,  that  at  frequencies  smaller  than  1-2  kc/s,  these  are  barely 
precise,  since  both  the  leakage  of  these  waves  through  the  ionosphere, 
and  the  spectrum  of  the  source  itself  at  f  <  i_2  kc/s,  should  be  taken 
into  account.  (See  Supplement  7. ) 

Figure  264  gives  the  signals,  E(t,  r),  computed  in  this  manner. 
Their  upper  part  presents  photo-oscillograms  of  atmospherics 
which  were  most  frequently  received  at  these  distances.  An  ex¬ 
ceptionally  good  correlation  between  the  oscillograms  and  the  results 
of  calculation  is  plain,  tending  to  establish  the  fact  that  the  theory 
apparently  takes  into  account  the  basic  peculiarities  of  propagation 
of  the  spectrum  of  the  waves  which  form  atmospherics. 


t 

Eqn.  (31.  7)  is  only  valid  if  z  k^  ^^<<  1,  (JRW). 
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Fig.  262  -  Product  of  moduli  of  the  interference  multiplier  and 
the  spectral  density  of  the  atmospheric  at  various  distances  (in  km). 


WO 


I  to  W'  10' 

Frequency,  c/s 

Fig.  263  -  The  sum  of  supplemental  phase  and  of  the  argument  of  the 
spectral  density  of  the  atmospheric  various  distances  (in  km). 
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Fig.  264  •  Comparison  of  observed  atmospherics  with  the 
computed  theoretical  signals. 


Experimental  comparison  of  the  modulus  of  the  spectral  density 
of  atmospherics  with  theoretical  computations  also  gives  a  good 
correlation  (Figs.  265  and  266).  However,  the  most  important  re¬ 
sult  of  harmonic  analysis  of  atmospherics  is  the  determination  of 
average  phase  velocity  of  the  electromagnetic  waves  over  the  range 
of  low  frequencies  (Fig.  267).  The  most  frequently  encountered 
experimental  values  of  v/c  ,  given  in  Fig.  267a,  coincide  well  with 
the  theoretical  values  v/c  up  to  where  frequencies  are  between  2  and 
3  kc/s.  At  lower  frequencies,  the  theoretical  value  of  velocity 
gradually  decreases  and  becomes  smaller  than  c,  but  experiments 
show  a  growth  in  velocity.  As  mentioned  above,  this  is  apparently 
explained  by  the  fact  that  leakage  through  the  ionosphere  becomes 
important  at  very  low  frequencies.  In  theory,  this  is  not  taken 
into  account. 


o 


Fig.  267-Theoretical  (1)  and 

(v/c)  a 
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Fig.  26 5 -Experimental  (a)  and  theo¬ 
retical  (b)  dependence  of  relative 
value  of  modulus  of  spectral  density 
of  atmospheric  on  various  distances 
(km)  in  England. 


Fig.  266 -Experimental  (a)  and  theo¬ 
retical  (b)  dependence  of  relative 
^  value  of  modulus  of  spectral  density 
of  atmospheric  on  various  distances 
(km).  Measurement  was  conducted 
in  Moscow. 


On  level 

0.  8 

0.  5 

Most 

probable 

values 


-43- 


Using  experimental  values  of  v/c  permits  us  to  determine,  with 
the  aid  of  the  theoi  etical  dependence  v/c  on  N/v,  the  values  of 
(N/v)g£f  observed  under  real  conditions.  In  the  given  experiments, 
the  dependence  of  N/v  on  frequency,  presented  in  Fig.  267b,  was  ob¬ 
tained.  Naturally,  the  most  probable  values  of  N/v  coincided  with 
those  given  in  computations  since  there  is  an  accord  between  corres¬ 
ponding  values  of  v/c  . 

It  should  be  noted  that  in  the  more  thorough  analysis  of  signals, 
E(t,  r),  which  were  computed  theoretically  (see  Fig.  264),  it  appears 
that  for  distances  of  2000-3000  km  and  greater,  the  beginning  of  the 
signal  is  somewhat  different  than  in  the  observed  atmospherics.  It 
is  possible  that  this  is  linked  to  the  fact  that,  in  the  computation  of 
£(t,  r),  the  finite  conductivity  of  the  earth  and  its  sphericity  were  not 
taken  into  account.  For  the  spectrum  of  waves  under  consideration, 
this  would  be  important  at  greater  distances  and  higher  frequencies. 
The  higher  frequencies  are  responsible  for  the  beginning  part  of  the 
atmospheric.  In  conclusion,  it  should  be  noted  that  the  atmospherics 
investigated  here  form  only  one  part  of  these  signals,  which  may  be 
termed  the  high  frequency  part,  and  which  have  a  duration  of  the  order 
of  0.  5  -  1  msec.  However,  behind  this  part  of  the  signal  there  fre¬ 
quently  follows  a  "tail”  where  duration  reaches  10-15  msec,  and  more 
(Fig.  268).  This  part  of  the  signal  is  formed  mainly  by  waves  whose 
frequencies  are  less  than  500-1000  c/s.  At  present,  theoretical 
computations  of  the  form  of  signal  in  this  range  have  not  been  carried 
out  to  an  adequate  extent.  Also,  the  "tails"  of  atmospherics  have  not 
been  the  subject  of  much  experimental  study. 


Fig,  268  -  Oscillogram  of  an  atmospheric 
with  a  low  frequency  tail. 
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2.  Whistling  Atmospherics 

The  whistling  atmospheric  is  received  by  ear  as  a  signal  whose 
frequency  changes  with  time.  At  first,  a  low  tone  is  heard;  then  a 
high  tone,  t  Most  frequently,  signals  of  such  type  run  from  400-8000 
c/s  [174,  175],  Apparently  the  greatest  intensity  of  the  whistling 
atmospherics  is  in  the  frequency  range  3000-4000  c/s. 

Records  of  a  single  whistling  atmospheric  and  a  group  of  such 
signals  are  presented  in  Figs.  269  and  270.  The  analysis  of  such 
signals  has  shown  that  they  are  associated  with  lightning  discharges 
originating  at  points  of  earth  which  are  connected  with  the  point  of 
observation  by  a  (terrestrial)  magnetic  line  of  force  (Fig.  271).  This 
type  of  whistling  atmospheric  is  called  a  short  one.  It  is  late  in  re¬ 
lation  to  the  lightning  discharge  with  respect  to  its  time  of  propagation 
along  the  trajectory.  However,  atmospherics  axe  also  received  which 
are  associated  with  local  thunderstorms.  These  propagate  along  a 
double  path;  that  is,  they  come  back  to  the  point  of  observation  after 
being  reflected  from  the  conjugate  point  of  the  earth  (Fig.  271).  In 
this  case,  the  signal  is  called  a  long  one,  since  the  transition  over 
the  whole  range  of  frequencies  (that  is,  the  whistle),  is  accomplished 
within  a  double  time.  The  time  of  retardation  of  the  beginning  of 
whistling,  in  relation  to  the  lightning  discharge,  is  directly  determined 
when  the  long  signals  are  received.  By  this  method,  the  time  of  group 
retardation  of  these  signals  was  established  as  of  the  order  of  2-3  sec. 
Furthermore,  it  follows  that  their  trajectory  reaches  lengths  oi  many 
thousands  of  kilometers  (see  below). 


t 

Usually  the  high  tone  precedes  the  lower  tone  in  conventional 
whistlers  (JRW). 
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Storey,  L.  R.  ,  An  investigation  of  whistling  atmospherics,  Phil. 
Trans.  Roy.  Soc.  ,  246,  113  (1953);  Can  J.  Phys.  ,  M,  153  (1956); 
3£,  1107  (1957). 
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Iwai,  A.  ,  Outsn,  J.  ,  On  an  investigation  of  whistling  atmospherics 
in  Japan,  Proc.  Res.  Inst.  Atm.,  4,  29  (1956). 
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presentation  of  trajectory  of  the  short  and  long 
(whistling)  atmospheric. 
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Experimental  data  have  made  it  possible  to  establish  that  the 
change  of  tone  (frequency)  of  the  signal,  with  time,  is  well  described 
by  the  formula  _ 

n/T  =  -  .  (31.9) 


The  quantity,  D,  is  named  the  dispersion  of  the  whistling  atmospheric. 
In  various  cases,  for  short  atmospherics,  it  changes  from  approxi¬ 
mately  10  sec  to  100  sec  and  more.  Analysis  of  the  sum  of  the  experi¬ 
mental  data  has  shown  that  the  dispersion  of  long  signals  is  approxi¬ 
mately  twice  as  large  as  the  dispersion  of  short  signals,  and  that 
there  is  some  correlation  between  probability  of  appearance  of 
whistling  atmospherics,  the  state  of  the  atmosphere,  and  magnetic 
activity.  This  lets  us  suggest  that  this  observed  phenomenon  can  be 
explained  in  the  following  way  [174] : 


As  already  pointed  out  in  the  preceding  paragraph,  the  waves 
radiated  by  a  lightning  discharge  partially  leak  through  the  ionosphere, 
especially  at  frequencies  between  10^  -  10^  c/s.  From  an  analysis  of 
the  coefficient  of  refraction  of  the  ionosphere,  taking  into  account  the 
magnetic  field  of  the  earth,  it  follows  that  this  occurs  for  the  quasi¬ 
longitudinal  propagation  of  these  waves  in  the  ionosphere.  At  these 
frequencies,  when  angles  as  small  as  5-10°  obtain  between  the  direction 
of  propagation  of  the  wr.ve  with  a  vector  of  the  magnetic  field,  H^, 
the  coefficient  of  refraction  becomes  imaginary.  In  this  way,  the 
waves  passing  through  the  ionosphere  are  canalized  (ducted)  along  the 
magnetic  field  of  the  earth,  and,  with  the  exclusion  of  the  case  where 
propagation  is  precisely  longitudinal,  are  regular  waves.  (This 
special  case  requires  a  detailed  supplemental  analysis  for  waves  of 
low  frequency.  )  Given  the  condition  that  everywhere  along  the  tra¬ 
jectory 

4rTN2 

— >>  1  (31.10) 

m  w 


(neglecting  the  influence  of  number  of  collisions),  the  coefficient  of 
refraction  of  these  waves  is  given  by  (see  1.  8)): 
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(31.11) 
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Storey,  L.  R.  ,  An  investigation  of  whistling  atmospherics. 
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where  6  is  the  angle  between  the  normal  to  the  wave  front  and  the 
magnetic  field  vector  H.  Computing  the  time  of  group  retardation  of 
signals  ^ 

t'  =  \  —  dS  ,  (31.12) 


u 


where  S  is  the  trajectory  of  the  wave,  a  is  an  angle  between  the  group 
velocity  vector  and  tangent  to  the  trajectory;  then 


u 


dn 

n  +  0)  — 
d  u 


(31.13) 


[see  (1.28)],  determines  group  velocity. 

Utilizing  Lhe  group  velocity,  it  is  possible  to  obtain  for  the  case 
when  th-  longitudinal  component  of  gyroscopic  frequency  equals 

G  H 

-  cos  B  »  uj  ,  (31. 14) 

me 


*k  (31.15) 

(s) 

Comparing  (31. 15)  with  (31.9),  we  see  that  the  expression  in  the  braces 
is  the  quantity  of  dispersion,  D.  If  the  condition  (31.  14)  is  not  ful¬ 
filled  everywhere  along  the  trajectory  of  the  wave  {jacket,  then, 

'  T  JiT^e  [-"ZnfrTc]  ° 

^  (8)  1  - 

and  it  is  possible  to  show,  in  distinction  to  the  first  case,  when  it  is 
possible  to  utilize  formula  (31.15),  that  the  time  of  retardation  of  the 
signal  does  not  increase  constantly  with  decrease  of  frequency.  There 
is  an  extreme;  at  some  value  of  frequency,  the  time  of  retardation 
reaches  a  minimum.  On  both  sides  of  this  value,  f,  time,  t,  in¬ 
creases  for  waves  of  greater  frequency,  and  also  for  waves  of  lesser 
frequency.  It  is  difficult  to  detect  this  effect  experimentally,  because 
under  usual  conditions,  the  border  value  of  frequency  lies  in  the 
region  of  waves  which  seep  through  the  ionosphere  badly.  However, 
here  are  some  experimental  proofs  of  this  effect;  the  extreme  may 
be  observed  in  some  series  of  experiments. 
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The  theoretical  explanation  of  whistling  atmospherics  given 
above  cannot  pretend  to  be  a  full  description  of  this  phenomenon. 

To  date,  there  is  no  sufficiently  accurate  method  of  taking  into 
account  the  influence  of  ions.  This  effect  may  be  substantial  at  low¬ 
er  frequencies  in  this  range.  Also,  the  computations  for  the  traject¬ 
ory  were  not  made  with  necessary  completeness  and  precision,  nor 
account  taken  of  anisotropy  of  the  ionosphere,  and  so  on.  .  . .  There¬ 
fore,  many  important  properties  of  this  phenomenon  are  not  made 
clear.  One  paper  (Fig.  272)  [176],  indicates  that  the  trajectory  of 
whistling  atmospherics  should  not  coincide  with  the  lines  of  the 
(terrestrial)  magnetic  field.  Therefore,  it  follows  that  the  point 
linked  to  the  source  does  not  lie  at  the  end  of  the  corresponding  field 
line.  However,  in  spite  of  the  absence  of  theoretical  explanation  of 
this  phenomenon,  it  is  still  possible  to  obtain,  with  the  aid  of  the 
given  approximate  formulae,  some  information  on  density  of  the  medi¬ 
um,  path  along  the  magnetic  field  of  the  earth,  and  other  quantities 
along  the  trajectory  of  whistling  atmospherics.  These  regions  stretch 
to  distances  as  remote  as  10,000-15,  000  km  from  the  earth  and  more 
(see  Fig.  272). 

In  conclusion,  it  is  expedient  to  note  the  following:  An  attempt 
was  made  to  ascertain  whether  effects  similar  to  the  whistling  atmo¬ 
spherics  can  be  observed  through  radiation  of  radio  waves.  Recently, 
the  results  of  observations  of  similar  echo-signals  at  a  frequency  of 
15.  5  kc/s  [177]  were  described  for  the  first  time.  These  signals  were 
retarded  for  approximately  1  sec.  relative  to  the  moment  of  their 
radiation.  From  analysis  of  experimental  data,  the  authors  of  the 
quoted  paper  came  to  the  conclusion  that  they  were  watching  signals 
propagating  along  the  line  of  force  of  the  magnetic  field  of  the  earth. 
Direct  proof  of  the  effect  of  whistling  atmospherics  with  the  aid  of 
radio  waves  was  thus  obtained  for  the  first  time.  It  is  possible  to 
suppose  that  this  phenomenon  played  a  part  in  various  earlier  radio 
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Maeda,  K.  ,  and  I.  A,  Kimura,  A  theoretical  investigation  on  the 
propagation  path  of  the  whistling  atmospherics.  Rep.  Ion.  Res 
Jap.  ,  ^,  No.  3,  10  5  (1956). 
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Helliwell,  R.  A.  ,  and  E.  Gehrels,  Observations  of  magneto -ionic 
duct  propagation  using  man-made  signals  of  very  low  frequencv 
Proc.  IRE,  No.  4,  785  (1958). 
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observations;  however,  proper  attention  was  not  given  it. 


Fig.  272-  Trajectories  of  propagation  of  whistling 
atmospherics  (solid  lines);  magnetic  force  lines 
are  presented  by  dots. 
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Appendix  6 

The  Influence  of  Sphericity  of  Earth  on  the  Problem  of  Wave -Directed 
Propagation  of  Electromagnetic  Waves  of  Low  Frequency 

1.  In  paper  [169]  ,  the  following  approximate  equation^  of  poles  was 
obtained  for  the  homogeneous  spheric  earth,  characterized  by  com¬ 
plex  coefficient  of  refraction  e*(ki  =kQ  and  homogeneous  con¬ 

centric  ionosphere,  characterized  by  the  coefficient  of  refraction, 

P2(C): 

k® 

i  ^  -  k^  =  fi  k  C  +  (2  R  )''] 

kf  ^  o  L  o  n  o  J 

[l-Pi(C  )  exp(i  2k  z  C  )][1  +  i(S  /C  )2(2  R  k  C  )'^  ]  .  (D.  28) 
n  on  nn  oon 

^  ,  -1 
[1  +  p-(C  )  exp(i  2k  z  C  )]  [1  -  i  (S  /C  )®(2R  k  C  )  ] 
n  on  nn  oon 

Here,  R  is  the  radius  of  earth;  z  is  altitude  of  the  ionosphere; 

Cn  are  poles;  =  \^  -  and  k^  =  cj/c.  This  equation  is  presented 
for  approximate  boundary  conditions 

^  [R  MR)]  =  -  i  ^  M  (R  y  R))  (D.  29) 

when  R  and  f^  are  radial  functions,  satisfying  the  equations 

d#  <R '„)  ^  ['''(R)  - 

It  is  easy  to  be  convinced  that  at  R  —  oo  (flat  earth)  the  equation 
(D.  28)  passes  into  the  equation  of  poles  for  the  flat  earth  (30.4). 
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Bremmer,  H.  National  Bureau  of  Standards  Rep.  No.  5518  (1957). 

t 

It  should  be  noted  that  Al'  pert  here  uses  an  exp  (  -  i  w  t)  whereas, 
in  (D.  31)  and  elsewhere  in  the  book,  exp  (i  w  t)  is  employed  (JRW). 
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The  equation  (D.  28)  is  rather  complicated  by  itself  and  is  not 
solved  in  [169].  However,  with  its  help,  and  can  be  expanded 
into  Taylor' s  series  by  powers  of  1/R  in  the  vicinity  of  1/R  -  0,  to  de¬ 
termine  the  influence  of  the  sphericity  of  the  earth.  Then  we  adopt 
supplemental  components,  A  S^,  A  Cj^,  to  the  values  of  Sjj  and 
for  the  flat  earth,  by  computation  of  the  corresponding  derivatives 

8C  88  C 

_ n  _ n 

»<4>  ’ 

and  so  on,  of  the  right  part  (D.  28).  As  a  result,  one  obtains  a  rapid¬ 
ly  convergent  series  in  which  it  is  sufficient  to  take  into  account  only 
the  first  term  of  the  series.  Similar  computations,  given  in  [243] , 
have  shown  that  in  many  cases  the  values,  A  and  A  C^,  are  quite 
large,  commensurate  with  and  even  larger  than,  and  C,^.  Apparent¬ 
ly,  it  is  possible,  using  this  method,  to  compute  sphericity  of  the 
earth,  accurately  based  in  each  separate  case  on  those  conditions 
suitable  for  the  formula  (D.  28),  and  using  the  series  for  and  Cj^. 


2.  In  a  recently  published  article  [242],  a  substantial  account  is 
given  of  the  present  state  of  the  theory  of  propagation  of  electro¬ 
magnetic  waves  of  low  frequency.  This  work  also  gives  the  results 
of  a  corresponding  analysis  of  the  spherical  waveguiding  problem. 
When  the  condition^  , 

c|  <<  1  ,  (D.  30) 
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Bremmer,  H.  ,  National  Bureau  of  Standards  Rep.  No.  5518  (1957). 

243 

Borodina,  S,  V.  ,  Yu.  K.  Kalinin,  G.  A,  Mikhailova  and  D.  S.  Fligel,  A 
survey  of  the  present  state  of  investigations  of  the  propagation  of  very 
long  electromagnetic  waves.  Radiophysics  (Bui.  Higher  Educational 
Institutions)  3,  6  (I960). 

242 

Wait,  J.  R.  ,  Terrestrial  propagation  of  very  low-frequency  radio 
waves,  a  theoretical  investigation,  J.  Research  NBS  64D  (Radio 
Prop.),  No.  2,  153-204  (Mar. -Apr.  I960). 

t 

In  eqn.  (D.  30)  <  <  should  be  >>.  This  was  a  misprint  in  reference 
[242].  For  frequencies  in  the  range  8-30  kc/s  this  condition  is 
usually  violated  for  the  dominant  modes  (JRW). 
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is  fulfilled,  the  electromagnetic  field  can  be  computed  with  the  aid 
of  formula  (30. 1),  or  (30.6),  with  the  multiplier,  0/sin  d  ,  of  geo¬ 
metric  character,  which  takes  into  account  the  correction  for 
sphericity  (see  above).  Here,  however,  the  poles,  C^,  and  corres¬ 
pondingly,  wave  numbers,  S^,  are  determined  from  an  equation 
of  poles  of  the  form 

^  z  1 

p,(C__)  p,(c;_)  expj.  i  r[c»  +  f  sj  ‘dx|  =  e-  ,  (D.31) 

o  ° 

where,  supplementing  designations  used  in  par.  30 

Ck  =  -  S''  ■  -  S'  =  • 

O 

The  coefficients  of  reflection  from  the  earth  and  the  ionosphere 
are  described  as  before  with  the  aid  of  Fresnel*  s  formulae.  When 
computing  the  general  solution  of  the  spherical  problem,  with  a  high¬ 
er  degree  of  precision,  when  limitation  (D.  30)  is  not  fulfilled,  the 
equation  of  poles  has  a  considerably  more  complicated  form.  This 
is  not  treated  here;  refer  to  paper  [242],  Also,  there  is  no  extended 
comment  on  the  results  of  analysis  of  the  corresponding  solution  of 
the  waveguiding  problem  for  the  region  close  to  the  antipode  [9,  163, 
242];  nor  is  there  an  examination  of  the  effect  of  the  magnetic  field 
of  the  earth  [167,  242];  and  many  other  questions. 


242 

Wait,  J.R.  ,  Terrestrial  propagation  of  very  low-frequency  radio 
waves,  a  theoretical  investigation. 

9 

Al*  pert,  Ya.  L.  ,  On  the  propagation  of  low  frequency  electromagnetic 
waves  over  the  earth*  s  surface,  Acad.  Science,  USSR (1954);Progre8s 
in  Physics,  60,  No.  360  (1956) 

163 

Schumann,  W.  ,  Uber  die  Oberfelder  bei  der  Ausbreitung  langer 
elektrische  Wellen  in  System  Erde-Luft-Ionosphare. 

167 

Budden,  K.  G.  ,  The  propagation  of  a  radio  atmospheric. 
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Appendix  7 

Passage  of  Electromagnetic  Waves  of  Low  Frequency 
Through  the  Ionosphere 

For  the  purpose  of  computing  the  influence  of  leakage  of  electro¬ 
magnetic  waves  through  the  ionosphere  (upper  wall  of  the  waveguide), 
it  is  necessary  to  compute  more  accurately  the  coefficient  of  re¬ 
flection  of  electromagnetic  waves  of  low  frequency  from  the  iono¬ 
sphere,  Pg  (C),  taking  into  account  the  influence  of  the  magnetic  field 
of  earth.  This  is  very  important  for  the  lower  frequencies  under 
consideration,  and  is  necessary  to  perform  for  both  the  spherical 
and  the  plane  problem.  Such  computation,  which  is  somewhat  com¬ 
plicated,  is  not  realized  practically  in  any  of  the  quoted  works.  How¬ 
ever,  it  is  expedient  to  point  out  here,  there  is  a  possibility  of  making 
a  more  simple  computation  of  poles,  based  on  the  physical  substance 
of  the  corresponding  solution  of  this  problem,  presented  in  Sec.  30. 

For  this  purpose,  it  is  necessary  to  solve  the  equation  of  poles, 
introducing  instead  of  the  coefficient  of  reflection,  pg(C),  the 
quantity  <» 

Pg(C)  [1  -  exp  i  k  ^  (n  -  i  x)  dz]  ,  (D.  32) 

O  V 

o 

where  ^ 

exp  i  k  \  (n  -  i  x)  dz 

o 

is  a  function  characterizing  the  transmission  into  the  ionosphere  of 
waves  of  the  corrresponding  frequency. 


Since  only  the  extra-ordinary  wave  in  quasilongitudinal  propaga¬ 
tion  passes  through  the  ionosphere,  it  is  possible  to  use  the  co¬ 
efficients  of  refraction  and  attenuation,  n,  and,  x,  given  above  in 
Supplement  1  [see  (D.  6)],  for  compu^ation  of  (D.  32).  The  results  of 
computations  of  the  transmission  function  of  the  ionosphere  for  a  real 
model  of  the  upper  atmosphere  (s3e  Sec.  6)  are  given  in  Table  D; 
sep>arately  for  day  and  night.  Corresponding  numerical  integration  is 
conducted  to  the  point  where  z  =3000  km  (upper  limit  of  integration) 
which  is  sufficient  for  these  evaluations.  In  Table  4  values  are  given 

3000 


exp 


X  dz 
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and  the  relation  v/c, of  average  phase  velocity  of  wave  v  to  c, 
characterizing  modulus  and  argument  of  function  for  passage  through 
the  ionosphere  by  ordinary  (index  +)  and  extra-ordinary  (index  -) 
waves.  At  present,  computations  of  poles  and  field,  taking  into  con¬ 
sideration  the  transmission  function,  are  not  yet  concluded.  There¬ 
fore,  there  is  no  possibility  of  comparison  with  experimental  results. 
Some  preliminary  evaluations,  however,  using  this  approach,  have 
shown  close  agreement  with  experimental  results. 


Table  4 


Transmission  function  of  the  atmosphere. 


Frequency 

Day  1 

1  Night 

(^)- 

i^l 

1 

i 

0,. 

1  9-10-* 

i  1,8-10-* 

1 

0,62  I 

2,8-10-3 

10 

3,8-10-* 

1,5-10-3 

;  2,1-10-3 

0,56  ! 

3,6-10-3 

50 

0,lu 

i  .>.io-">* 

3,10-3 

0,72 

4,8-10-3 

10* 

0,18 

~10-'«* 

4,10-3 

!  0,79 

5,3-10-3 

10» 

0,26 

_10-2.6-10* 

it. 8- 10-3 

0,83 

1,6-10-* 

2. 10* 

0,16 

__ 

— 

— 

j.lO* 

7,3-lU-* 

— 

~ 

— 

— 

7-10* 

8,3-10-* 

— 

— 

— 

— 

10* 

2,9-10-* 

_io-«-io* 

3,4-10-2 

1  0,59 

5,0-10-* 

105 

1,2-10-5 

;  l,0-10-» 

0,18 

0,14 

10“ 

10-'"' 

-^lO-'l'** 

0,17 

—  lo-'’-  )o< 

'\5o 

i 


It  should  be  noted  here  that  at  frequency,  f  <  10  c/s,  the  average 
phase  velocity  is  smaller  than  velocity,  c,  by  two  magnitudes  and 
more,  a  fact  which  is  linked  with  the  magnetohydrodynamic  character 
of  waves  at  these  frequencies. 
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Final  Notes 


As  a  matter  of  interest  to  readers  of  this  translation,  a  review 
of  an  earlier  paper  of  Al'pert  is  included  here.  This  was  reprinted 
from  Math.  Revs.  ,  vol.  21,  No.  10,  1285  (Nov.  I960).  Also,  I  am 
attaching  copies  of  letters  written  to  Borodina  and  Al'pert.  These 
have  never  been  answered. 

Al'pert,  Ya.  L.  Computation  of  the  field  of  LF  and  VLF 
radiowaves  over  the  Ei^h’s  surface  under  natural  condi¬ 
tions.  Radiotehn.  i  Electron.  1  (1956),  281-292;  avail¬ 
able  in  translation  T320  by  E.  R.  Hope.  Directorate  of 
Scientific  Information  Service,  Defense  Research  Board, 

Ottawa,  Canada,  Nov.  1959. 

The  paper  sets  out  to  compute  the  time-harmonic 
electromagnetic  fields  in  the  frequency  range  500  c/s  to 
.*10  kc/s  at  large  distances  along  the  surface  of  the  earth, 
from  a  transmitting  antenna.  The  adapted  model  is  a 
))arallel-plate  waveguide  whose  lower  boundary  is  the  air- 
ground  interface  and  whose  upper  boundary,  at  height  h, 
is  the  lower  edge  of  the  ionosphere  D  region.  'I'he  source  is 
represerded  by  a  vertical  electric  dipole.  The  calculations 
arc  based  on  the  general  integral  contour  representation  - 
fur  the  ])lanar  )>rublem.  'fhis  is  transformed  to  a  sum  of 
inodes  u’hich  individually  correspond  to  the  residues  of  the 
jioles  of  the  integrand.  'Flie  resulting  pole-determining 
e(|uation  involves  the  reflection  coefficient  />2(C)  of  the 
upper  boundary  for  ]>lane  waves  incident  at  a  (complex) 
angle  whose  cosine  isC.  'I'he  conductivity  of  the  D  region 
is  allowed  to  vary  with  height  h  +  z  according  to  the 
function  [  I +e.\p  ( -/»:)]  '  where  m  is  chosen  to  be 
to  corn'spond  to  be.st  available  jmblished  data  (the 
iiiHuenee  of  the  earth's  magnetic  field  is  neglected). 

Taking  1*.  S.  Epstein's  formula  [Proc.  Nat.  Acad.  Sci. 
l’.S..A.  16  (1930),  627-037)  for  the  reflection  coefficient 
the  author  then  u.ses  a  .seni’-empirical  technique 
described  as  tiie  “alignment  methocr’,  to  estimate  g2(C). 

This  is  the  crucial  step  in  the  analysis  and  its  implications 
are  not  a<le<juately  diseu.ssed  in  the  imticr.  Actually, 

Epstein's  expres.siou  for  g(t')  is  valid  only  for  horizontal 
polarization  (i.e.,  the  electric  vector  is  horizontal),  whereas 
the  reejuired  e.xpre.ssiou  for  the  pole  eipiation  corre¬ 

sponds  to  vertical  iiolarization  (i.e..  the  magnetic  vector  is 
horizontal).  'I’he  “alignment  "  procedure  amounts  to  using 
Epstein's  ivsults  for  horizontal  jiolarization  to  define  an 
equivalent  sharply  boundeil  ionosphere  whose  conduc¬ 
tivity  is.  strictly  s|»eaklng.  both  a  function  of  C  and  the 
frequency.  It  is  then  a.ssunied  that  this  same  “effective'’ 
conductivity  may  be  inserted  into  the  usual  Fresnel 
reflection  formula  for  vertical  polarization  to  yield  a 
numerical  value  for  pi(C).  In  view  of  the  radically  different 
behavior  of  the  Fresnel  reflection  coefficients  for  horizontal 
and  vertical  ])olarizntion  at  grazing  angles  (i.e..  small  C) 
the  validity  of  this  proce<lHre  is  doubtful. 

'fhe  author  then  claims  that  his  results  may  be  corrected 
for  the  s|)hericity  of  the  earth  by  multiplying  the  fields  by 
the  factor  (fl/sin  0)'  -  where  6  is  the  central  angle.  The 
adequacy  of  this  stejj  is  questioned  since  the  earth 
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ciirvatiirc  inodiHcs  tlic  {Mjlc-deteriniiiiiig  equation  and  the 
attiMuiation  of  the  dominant  mcaics  would  be  increased  if 
this  were  taken  into  account.  This  fact  would  have  been 
evident  if  the  ])ioblem  had  Ir'cii  formulated  in  sjjherical 
•ieometry  at  the  outset.  Actually,  an  earth-flattening 
apj»ro.\imation  for  the  )H)le-determining  equation  is  only 
valid  when  the  spherical  wave  functions  are  adequately 
reju’esented  by  t  heir  second  order  or  Debye  appro.ximation. 
Such  is  the  ca.se  when  1,  where  k  is  the  wave 

number  aiul  «  is  the  radius  of  the  earth.  This  ap[)ears  to  be 
violated  for  many  of  the  numerical  results  given  in  the 
])aper.  In  fact,  if  curvature  is  completely  neglected,  a 
further  restriction  (which  couhl  l>e  obtained  from  geo¬ 
metrical  optics)  is  that  1.  'I'his  is  also 

violated  in  many  of  the  ea.ses  treated  by  the  author. 

It  is  possible  that  the  re|M)rted  agreement  with  e.xperi- 
mcntal  data,  in  the  frequenc\'  range  Jo  to  .‘{0  kc/s,  is 
jiartly  fort\utous  since  the  "alignment  j)rocedure"’, 
mentioned  above,  yields  overly  high  attenuation  of  the 
nodes,  whereas  the  neglect  of  earth  curvature  leads  to  very 
ow  attenuation.  J.  J{.  Waif  (Boulder,  Colo.) 
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U.  S.  DEPARTMENT  OF  COMMERCE 
NATIONAL  BUREAU  OF  STANDARDS 
BOULDER  LABORATORIES 
BOULDER.  COLORADO 

March  3,  1961 


CRYOOENtC  ENQINCERING  LABORATORY 
CENTRAL  RADIO  PROPAGATION  LABORATORY 
RADIO  standards  LABORATORV 

AIR  MAIL 


(N  RBRLV 

REFER  TO  FILE  NO. 


79.  70 


#Dr.  S.  V.  Borodina 

Institute  of  Terrestrial  Magnetism,  Ionosphere  and  Propagation 
Moskovskaya  obi.  ,  Leninskiu  r-n 
P/O  Vatutenki,  USSR 

Dear  Dr.  Borodina: 

I  have  been  very  interested  in  your  paper  "Review  of  the  Modern 
State  on  the  Propagation  of  Ultralong  Electromagnetic  Waves,  "  co¬ 
authored  by  lu.  K.  Kalinin,  G.  A.  Mikhailova,  and  D.  S.  Fligel.  I  have 
also  been  very  interested  in  the  papers  of  your  distinguished  colleague, 
Ya.  L.  Al'pert.  Our  laboratory  would  be  very  pleased  to  receive  re¬ 
prints  of  any  of  this  work.  We  would  also  like  a  copy  of  your  paper 
published  in  the  works  of  NIZMIR  13,  3  (1959). 

In  exchange  for  the  above  material,  I  am  sending  to  you  under 
separate  cover  twenty  reprints. 

I  was  rather  surprised  about  your  comments  on  mode  numbering. 
Isn't  it  really  a  matter  of  nomenclature?  It  is  immaterial  whether  the 
dominant  mode  is  labelled  zero  or  one,  provided  that  the  excitation  is 
properly  accounted  for.  Similar  remarks  apply  to  the  higher  mode. 
Some  further  discussion  on  this  point  by  Dr,  H.  H.  Howe  is  contained 
in  Vol.  IV  of  the  Proceedings  of  the  19  57  Symposium  on  VLF  radio 
waves.  A  copy  of  this  is  also  being  sent  under  separate  cover. 

I  believe  the  question  of  sphericity  of  the  earth  is  a  very  important 
factor.  My  view  on  this  is  expounded  in  some  of  the  reprints  being 
sent  to  you.  In  this  connection  I  think  you  might  also  be  interested 
in  my  review  of  Al'pert' s  article.  A  copy  of  this,  taken  from  Mathe¬ 
matical  Reviews,  is  enclosed.  I  would  appreciate  it  if  you  could  bring 
this  to  his  attention.  (I  wasn't  sure  of  his  current  address.  )  As  you 
see,  1  am  particularly  concerned  about  the  applicability  of  the  flat- 
earth  model  to  frequencies  above  15  kc/s  since  the  condition 
(ka/2)^' ^  C  >  >  1  is  seriousl/  violated.  I  would  be  very  interested 
to  know  if  you  agree  with  my  view  on  the  subject. 

Yours  truly, 


End  James  R.  Wait,  Consultant 

Rev  of  Al'  pert' s  Boulder  l  aboratories 

article 
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U.  S.  DEPARTMENT  OF  COMMERCE 
NATIONAL  BUREAU  OF  STANDARDS 
BOULDER  LABORATORIES 
BOULDER.  COLORADO 

June  15,  1961 

CRVOQENIC  ENOINEERING  LABORATORY 
CENTRAL  RADIO  PROPAGATION  LABORATORY 
RADIO  STANDARDS  LABORATORY 

AIR  MAIL 
* 

Dr.  S.  V.  Borodina 
Institute  of  Terrestrial  Magnetism, 

Ionosphere  and  Propagation 
Moskovskaya  obi. 

Leninskiy  r-n 
P/D  Vatutenki 
USSR 


Dear  Dr.  Borodina: 


I  was  wondering  if  you  received  my  earlier  letter  (dated  March  3, 
1961)  and  a  number  of  reprints.  A  copy  of  the  letter  is  enclosed.  I 
was  hoping  to  receive  some  comments  from  you  and  reprints  of  some 
of  your  papers. 

I  have  looked  again  into  the  question  of  the  influence  of  ground  con¬ 
ductivity  on  the  VLF  modes.  I  have  been  able  to  confirm  our  earlier 
calculations  which  were  quoted  in  Fig.  4  of  your  paper  in  Izvestia 
V.  U.  Z.  ,  Vol.  3,  pg.  12,  Jan.^eb.  I960.  This  was  done  by  carrying 
out  a  simple  perturbation  calculation  using  a  method  described  in  my 
paper  in  Proceedings  of  the  Inst,  of  Radio  Engineers,  vol.  45,  pg.  763, 
June  1957.  This  shows  that  for  low-order  modes  the  additional  attenu¬ 
ation  resulting  from  finite  ground  conductivity  is  given  by 


i 

6  14  2  3 

^  (G)  X  10^ 


1. 

decibels  per  1000  km  of  path  length,  where  G  =(e^  o/c)*.  For  example, 
if  h  =  80  km,  f  =16  kc/s,  a  =2  x  10”  mhos/m,  we  find  that  the  additional 
attenuation  due  to  ground  losses  is  1.6  db.^000  km.  This  agrees  very 
closely  with  the  difference  between  the  two  curves  for  c  =2.  08  and 
(T  =  00  in  Fig.  5  of  my  paper  on  pg.  770  of  the  same  issue  of  the 
journal.  These  latter  results  were  obtained  from  a  full  numerical 
treatment  of  the  exact  mode  equation  on  a  digital  computer. 


Dr,  S.  V.  Borodina 
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(2) 


June  15,  1961 


Could  there  be  an  error  of  amount  s/2  in  your  equation  (10)? 

This  might  resolve  the  discrepancy.  It  is  also  possible  that  in 
transferring  my  results  in  MKS  units  to  your  CGSE  system  an  error 
crept  in. 

I  hope  to  hear  from  you  shortly. 

Yours  truly, 


End  James  R.  Wait 

cy  Itr  3/3/61  Consultant 

Central  Radio  Propagation  Laboratory 
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U.  S.  DEPARTMENT  OF  COMMERCE 
NATIONAL  BUREAU  OF  STANDARDS 
BOULDER  LABORATORIES 
BOULDER.  COLORADO 

September  25,  1961 

CHYOOENIC  ENOINEBRINO  LABOEATOBY  im 

CENTEAL  EADIO  PEOEAOATION  LABOEATOEY  EEFEE  TO  FILE  NO. 

EAOlO  STANOAEDS  LABOEATOEY 

AIR  MAIL  79.  70 

* 

Dr.  Ya.  L.  Al'  pert 

Institute  of  Terrestrial  Magnetism,  Ionosphere 
and  Propagation 
Moskovskaya  obi. 

Leninskiy  r-n 
P/O  Vatutenki 
U.S.S.R. 


Dear  Dr.  Al'pert: 


I  have  been  reading  some  of  the  recent  Russian  literature  on  VLF 
propagation.  In  two  places^  you  and  your  colleagues  have  used  equation 
(25)  in  Bremmer' s  NBS  Report  5518  (later  published  in  NBS  Jour.  Res. 
Sec.  D.  ,  July  1959).  From  this  an  earth -curvature  correction  to  the 
flat -earth  mode  equation  has  been  obtained.  In  the  conventional 
notation  this  equation  would  read 


~  N/kf  -  k  S  =  li  k  C  + 

kf  ^  o  Lon  Cjj  2aJ 


[1  -  R.  (C  )  exp  (i  2k  h  C  )][1  +  i(S  /C  )2(2k  a  C  )  ] 
in  o  n  nn  o  n 

- =T" 

ri  +  R.  (C  )  exp(i  2k  h  C  )][1  -  i(S  /C  )*  (2k  a  C  )  ] 
'■  1  n  ^  o  n-*^  nn  o  n 


where  k^  and  k^  are  the  (complex)  wave  numbers  for  the  air  and  earth, 
respectively,  a  is  the  radius  of  the  earth,  his  the  height  of  the 
ionospheric  reflecting  layer  whose  reflection  coefficient  is  R£(Cj^), 
and  S^  are  the  usual  complex  and  cosine  and  sines  chaiacterizing 
the  mode.  Now,  I  am  wondering  if  this  equation  is  valid  since  the 
difference  between  t  and  t'  has  been  neglected.  Is  this  not  the  same 
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as  saying  that  C  and  C'  are  the  same?  Actually 
n  n 


Consequently,  should  not  2h/a  be  much  less  than  C®  to  allow  t  to  be 

replaced  by  t  '  ?  This  is  a  fairly  stringent  condition.  Another 

criticism  of  this  equation  is  that  the  indicated  curvature  correction 

can  only  be  significant  if  C®  k^a  is  not  large  compared  with  unity. 

But,  the  condition  for  the  use  of  the  asymptotic  forms  (eqns.  23  and  24) 

in  Dr.  Bremmer' s  paper  is  only  valid  if  k  a  is  large  compared 

with  unity.  Therefore,  I  question  whether  the  results  are  valid  under 

any  practical  condition.  Borodina  ,et  al.have  taken  this  equation  to 

obtain  curvature  corrections  to  the  attenuation  rates.  In  some  cases 

their  corrected  value  is  ten  times  larger  than  the  case  for  a  =  oo. 

Obviously,  this  results  from  the  smallness  of  C  ^  k  a. 

no 

I  would  be  interested  in  your  reply  to  my  questions.  Possibly, 
there  is  something  which  I  have  overlooked. 


Reprints  of  a  few  recent  papers  of  ovr  (e.  g.  ,  Wait  in  NBS  Jour. 
Res.  Sec.  D,  Mar. /A.pr.  I960  and  Jan. /Feb,  1961,  also  Wait  and  Spies 
in  J.G.  R.  ,  Aug.  I960  and  March  I960), are  enclosed.  In  return,  I 
would  appreciate  reprints  of  any  of  your  papers. 


Yours  truly. 


Ends  James  R.  Wait 

Reprints  143,145,161,169  Consultant 

Cy  Wait  review  of  Central  Radio  Propagation  Laboratory 

Al'pert  paper  in  Math 
Revs,  Nov.  i960 
Cy  Itr  dtd  6-15-61  to 
Dr.S.V.  Borodina 

P.  S.  Incidentally,  I  have  written  two  letters  to  Dr.  Borodina,  but 
as  yet  no  reply  has  been  forthcoming.  A  copy  of  the  last  one  is  en¬ 
closed.  Possibly  you  could  forward  this  to  her. 


